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ABSTRACT 

We investigate the X-ray properties of BzK-selected galaxies at z ~ 2 using deep X-ray data 
in the Chandra Deep Field South (CDFS) and Chandra Deep Field North (CDFN). A subset of 
these BzK galaxies have been proposed as Compton thick Active Galactic Nuclei (AGN) can- 
didates based on a high ratio of infrared (IR) and ultraviolet (UV) star formation rates (SFR). 
With the benefit of deep 24/im observations our sample of these IR excess galaxies is larger 
than previous studies, and combined with the deepest X-ray data yield new insights into the 
(-h \ population. We identify 701 and 534 star forming BzKs (sBZKs) in the range z = 1.2 — 3.0 in 

the CDFS and CDFN respectively. Of these we directly detect in X-rays 49 sBzKs in CDFS 
and 32 sBzKs in CDFN. Stacking the undetected sources also reveals a significant X-ray sig- 
nal. Investigating the X-ray detection rate and stacked flux versus the IR excess parameter (i.e. 
SFR tota ;/SFR;yy. corr ), we find no strong evidence for an increased X-ray detection rate, or a 
harder X-ray spectrum in IR Excess sBzKs. This is particularly the case when one accounts 
for the strong correlation between the IR excess parameter and the bolometric IR luminosity 
(Lm), e -g- when controlling for Ljr, the IR Non-Excess sBzKs show a detection rate at least 
as high. While both direct detections and stacking suggest that the AGN fraction in sBzK 
galaxies is high, there is no clear evidence for widespread Compton thick activity in either 
the sBzK population generally, or the IR Excess sBzK subsample. The very hard X-ray sig- 
nal obtained for the latter in earlier work was most likely contaminated by a few hard X-ray 
sources now directly detected in deeper X-ray data. The X-ray detection fraction of passive 
BzK galaxies in our sample is if anything higher than that of sBZKs, so there is no evidence 
for coeval black hole growth and star formation from X-ray analysis of the BzK populations. 
Because increased AGN activity in the IR excess population is not indicated by our X-ray 
analysis, it appears that the bulk of the IR Excess sBzK population are luminous star-forming 
galaxies whose SFRs are either overestimated at 24/im, underestimated in the UV, or both. 
This conclusion reinforces recent results from Herschel which show similar effects. 
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frared: galaxies - ultraviolet: galaxies - X-rays: galaxies 



1 INTRODUCTION 

Extensive deep multi-wavelength observations across the sky have 
vastly improved our knowledge of star-formation and AGN ac- 
tivity across a range of redshifts. It has become abundantly clear 
that the Universe at z ~ 1 — 3, was much more active than it 
is currently, at z ~ 0. For example, the star-formation rate den- 
sity steadily rises from z = to z = 1, flattens out u ntil at least 
z ~3 and then falling off again at the highes t redshifts dLillv et all 
ll996l:lMadau et alJl996l : lHopkins et alj2005h . AGN activity shows 
a similar form, indicating that globally, at least, the growth of super- 
massive black holes roughl y tracks the buildup of the stellar pop- 
ulation in galaxies (e. g . IBovle & Terlevicb] 1 19981 : Isilverman et al.l 
l2008l : lAirdetai]|2010h . Observing galaxies at z~2 should give us 
tremendous insight into this evolution, and the interplay between 
AGN and star-formation in their host galaxies during this more 
active phase, but spectroscopic identifications, and hence our un- 



derstanding of galaxies with 1.4 < z < 2.5 is complicated by 
the "redshift desert". At these redshifts there are very few spectral 
features at observed-frame optical wavelengths that can be used to 
determine the redshift. In response to this problem numerous se- 
lection techniques have been developed that use broad band pho- 
tometry instead of spectroscopy to select z ~2 candidate galax- 
ies. For example, the Lyman-break and Balmer bre ak features have 
been used to select ga laxies with 1 < z < 3 terbetalJliii 
lAdelberger etai]|2004h . 

The BzK selection technique of Dadd i et all J2004) combines 
optical (B and z band) and near-IR (K band) photometry to tar- 
get the Balmer and 4000 A break features of star-forming and pas- 
sively evolving galaxies at z ~2 (Reddy et al. 2005; Shapl ev et al.l 
l2005l : IWuvts et alj2008t) . The sele ction technique is expected to be 
largely independent of reddening ( Daddi et al. 2004) and the near- 
IR photometry allows extremely red passively evolving galaxies to 
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be selected, as well as bluer, star-forming ga laxies dThompson et al.l 
1 19991 ; iFranx et alj|2003l ; iMcCarfhvl 12004) . BzK selection shows 
good agreement with other 2 ~2 galaxy selection techniques, al- 
though it is liable to miss the bluest/younge st star-forming galax- 
ies (Kg > 21 Vega mag ; Redd v et alJ[20o3) . Many sBzKs show 
strong emission in the mid-IR at 2 4/xm, the majority o f which is 
attributable to dusty star formation Daddi et al. (2007 3). By com- 
paring U V and IR derived star-formation rates (SFR's) lpaddi et all 
d2007bh identified sBzKs with excess I R emission ( IR Exc ess 
sBzK) as candidates for containing AGN. Dad di et al.l d2007bl) to 
attempted to isolate sBzKs which might contain AGN using this 
mid-IR excess. AGN are expected to exhibit strong continuum 
emission at m id-IR wavelengths due to re-radiated emission from 
hot d ust (e.g. IPier & Krollkll 19921 ; lEfstathiou & Rowan-RobinsorJ 
1995), a feat ure that has been exploited previou sly to identify ob- 
scur ed AGN dLacv et al.ll2004IStern et alj2005h . 

iDaddi et alj d2007bl) stacked the X-ray emission of IR Ex- 
cess sBzKs in the GOODS-S field, yielding a hard X-ray sig- 
nal (r ~ 0.8^0 3) consistent with a significant population of 
Compton thick (Nh > 1.5 x 10 24 cm' 2 ) AGN. The stacked 
signal of X-ray undetected sBzKs that did not exhibit excessive 
IR emission (IR Non-Excess sBzKs) appeared softer (r ~ 1.7) 
and consistent with pure star-formation. These results have sig- 
nificant implications for co-eval star-formation and black hole 
growth and the composition of the X-ray background. The pres- 
ence of "buried" Compton thick AGN in star-forming galaxies 
(like sBzKs) h as been predicted by models of AGN-galaxy-co- 
evolution (e.g. Hopkins et alj 120051) . The majority of the X-ray 
background dGiacconi et alj 1 19621) has been res olved using the 
emiss ion of discrete X-ray sources (as predicted bv lSetti & Woltierl 
1989), but t here remains an u nexplained hard (E > 6 keV) X- 
ray deficit dWorslev et afll2005h . Models of the X-ray background 
predict the hard X-ray background deficit could represent a large 
popul ation of as yet undiscovered Compton thick AGN at 2 ~ 
1-2 dComastri et al.ll995l;lGilli. Comastri. & Hasingej2007l) . Ex- 
trapolating the stacked X-ray emission of the X-ray undetected 
IR-Excess sBzKs over all the sky would explain a substantial 
portion (~ 10% ~ 25%) of the hard X-r ay background deficit 
dDaddi et all l2007bl) . lAlexander et alj d201lh carried out further 
analysis of sBzKs using new, deeper 4Ms CDFS X-ray data and 
found that stacking X-ray spectra of heavily obscured (r < 1) 
sBzKs produced a hard, reflection dominated spectrum consis- 
tent with Compton thick AGN emission. The stacked spectrum 
features a large EW Iron line at rest-frame E ~ 6.4 keV simi- 
lar to that observed in Compton thick AGN in the local Universe 
dMatt. B randt. & Fabi anll996l;llwasawa & Comastrill9 98) and dis- 
tant Univers e ( Georgantopoulos et al. 2009; Comastri et alj 2011 ; 
iFeruglio et ail 1201 ll ; iGilli et alj 1201 ll) . lAlexander et alj J20T1I) 
stacked the spectra of sBzK galaxies containing heavily obscured 
AGN and found their spectra was best fit by a pure reflection model 
with a strong iron Kq line at rest-frame 6.4 keV. This is consistent 
with Compton thick absorption (Nh ~ 10 24 cm~2) in a small frac- 
tion of sBzK galaxies. Stacking X-ray undetected sBzKs in deeper 
4Ms CDFS data did not yield conclusive evidence of a high frac- 
tion of heavily obscured AGN in the more general sBzK popula- 
tion, however, some what in disagreement with the earlier work of 
Daddi et alj d2007bl) . These observations imply that the space den- 
sity of Compton thick AGN at z ~ 2 is comparable to that of 
unobscured AGN at 2 ~ 2. 

The IR excess analysis o f lDaddietal.l d2007b) is reliant upon 
the accuracy of IR and UV SFR estimates. Mid-IR observations 
have been found to be reliable proxies for bolometric IR lumi- 



nosities dDale et alj |2005|) and it has been common practice to 
fit 24/im observations (rest-frame 8/im at 2 ~ 2) with empir- 
ically derived star- formin g galaxy templates (e.g. ICharv & Elbazl 
l200ll ; |Pale & Heloull2002h to estimate the SFR. With the advent of 
Herschel our understanding of the Far-IR spectra of star-forming 
galaxies has impr oved markedly, esp ecially at z > 1. Using Far- 
IR Herschel data lElbaz et alj d201ll) calculated Lir for numer- 
ous galaxies with MIPS coverage and compared them to Lir 
estimates made using 24/im selected star-forming templates of 
ICharv & Elbazl d2001|) . Based on this analysis it is c oncluded that 
the star-forming templates of Charv & Elbazl j200Th are overesti- 
mating the bolo metric IR luminosity f or galaxies with Lir > 10 12 
Lq . Additionallv lNordon et alj d2012l) find that the IR Excess emis- 
sion at rest-frame 8/im can be explained entirely by PAH features, 
removing the necessity for mid-IR continuum emission from AGN 
assumed by previous IR selection techniques. The UV estimates of 
SFR are known to suffer in dusty (ie heavily star-forming) galax- 
ies because redden ing corrections becom e less accurate as the dust 
opacity increases (Gold aderet ail [2002). Thus, selection of AGN 
via comparison of IR and UV-derived SFRs may not be as ro- 
bust as previously thought. Indeed, recent wor k have suggested 
that some IR excess selection techniques (e.g. iFiore et aU 20081) 
have been shown to select equal parts heavily obscur ed AGN, low- 
luminosity AGN and purely star-forming galaxies (Donley et alj 
2008; Georgakakis et a l. 2010) when followed up using X-ray data. 

In the light of this, in the present paper we revisit the is- 
sue of the AGN co ntent in IR Excess BzK galaxies pioneered by 
lDaddietal.ld2007bl) . We examine the X-ray properties of IR lumi - 
nous sBzKs, selected in a similar manner to Da ddi et al . (2007b), 
using 4Ms CDFS and 2Ms CDFN X-ray data, through a mixture 
of direct detections and stacking. In Section |2~T| we discuss the op- 
tical and IR data used in this work. Section I2T2I outlines the BzK 
and 24/im selection methodology while Sections [23] and |2j4] out- 
line the X-ray data reduction and X-ray stacking technique respec- 
tively. In Section [3] we present the results of direct X-ray detec- 
tion of sBzKs and stacking of X-ray undete cted sBzKs and i n Sec- 
tion[4] we compare our findings to those o fPa ddietalj j2007b) and 
lAlexander et alj ( 1201 If ). In Section|5]we discuss the implications of 
these results and in Section|6]we summarise our findings. Through- 
out this work a standard, flat ACDM cosmology with Q.\ = 0.7 
and Ho = 70 km s~ x Mpc -1 is assumed. 



2 DATA 

2.1 Optical and IR data 

In the CDFS we use Hubble ACS F435W, F606W, F775W and 
F850LP band (henceforth B, v, i and 2 respectively), ISAAC K 
band, IRAC 3.6/iin and MIPS 24/xm band data from the GOODS- 
MUSIC-v2 catalogue dSantini et aU l2009h . All objects have ei- 
ther a spectroscopic (64%) or a photometric (36%) redshift. This 
dataset presents identical B, v, i, z, K and 3.6am photometry 
to the GOODS-MUSIC cat alogue dGrazian et ai]|2006l) that was 
used bv lDaddietal] d2007b) but incorporates deeper MIPS 24/im 
band data from the FIDEL survey (PI: Mark Dickinson). The 
24/im source fluxes in GOODS-MUSIC-v2 were identified and 
matched to the multi- w avelength data using the ConvPhot software 
dde Santis et alj |2007), as opposed to matching the MIPS sources 
to multi-wavelength data using a more basic fixed m atching ra- 
dius method (as was the case for IDaddi et"al] 120051 and subse- 
quently Paddi et al 2007a and 2007b). The additional depth and 
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Table 1. BzK galaxy samples and IR Excess and Non-Excess sBzK sub- 
samples for the CDFS and CDFN. 



Galaxy type 


Total CDFS 


Total CDFN 


Total Combined 


pBzK 


21 


5 


26 


nBzK 


2168 


2006 


4174 


sBzK 


851 


594 


1445 


IRX (1.2 < z < 3.0) 


158 


198 


356 


IRNX (1.2 < z < 3.0) 


543 


336 


879 



^9 




sBzK * 
pBzK X 
nBzK + 



2 3 
(B-z)„ 



Figure 1. BzK diagram of galaxies in the GOODS-MUSIC-v2 and MODS 
catalogues. Solid Black lines represent the boundaries between different 
BzK galaxy type criteria. 



use of the ConvPhot software allows the identification of additional 
24/im sources and t herefore increases our sample size compared to 
|Paddietal.N2007bl) . The GOODS-MUSIC-v2 catalogue has more 
spectroscopic redshifts than GOODS-MUSIC, improving the accu- 
racy of our luminosity calculations. In the CDFN we use Hubble 
ACS B, V, i and z band, MOIRCS K s , IRAC 3.6/im and MIPS 
24/im b and photometry collec ted for the MOIRCS Deep Survey 
(MODS: iKaiisawa et alj|201 ll) . MODS is the only publicly avail- 
able catalogue in CDFN to date to have matched K band to 24/im 
sources. Again, sources have a mixture of spectroscopic (70%) and 
photometric (30%) redshifts. 

To ensure the reliability of the BzK selection and IR excess 
analysis we require all the optically galaxies in our sample to satisfy 
the following criteria: 

(i) A secure K-band detection (Ks < 22 Vega). 

(ii) less than 0.5 arcseconds separation from an IRAC 3.6um 
selected source to verify that the source is genuine (large PSFs of 
24/im sources leads to source blending). 

(iii) Error on (B — z)ab < 0.4 to ensure reddening corrections 
are reliable. 

(iv) Either a spectroscopic or photometric redshift. 

In total 5,685 galaxies across CDFS and CDFN satisfy these criteria 
(plotted on BzK diagram; Figure[TJ. 
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Figure 2. Redshift distribution of all 1485 sBzK galaxies. Filled portion 
of the histogram represents sBzK galaxies subjected to IR Excess analysis 
(1.2 < z < 3.0), the remaining galaxies are omitted as they possess either 
z < 1.2 or z > 3. 
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Figure 3. Plot of SFRtotJSFRuv corr versus redshift for sBzK galaxies 
used for IR excess analysis. Galaxies above the horizontal dashed line sat- 
isfy the IR excess galaxy criterion, galaxies below are IR non-excess. Filled 
squares are galaxies with spectroscopic redshifts while open squares are 
galaxies with photometric redshift estimates only. Black crosses represent 
galaxies that have been directly detected in the 4Ms CDFS and 2Ms CDFN 
data. 



2.2 IR Excess sBzK selection 

BzK galaxies with excess IR emission compared to that expected 
for star formation, may have a significant contribution to the IR flux 
originating from accretion and are therefore candidate obscured 
A GN. To define IR Ex cess sBzK galaxies we follow the methods 
of iDaddi et at] d2007bl) . using a combination of the BzK selection 
technique foaddi et al.ll2004l) and SFR's calculated using UV and 
IR obs ervations. The sta r-forming BzK (sBzK) selection (Equation 
rn>of lDaddietal.N2004h identifies star-forming galaxies at z > 1.4. 



(z-K)ab - {B-z) AB > -0.2 



(1) 



© 0000 RAS, MNRAS 000, 000-000 



4 C. Rangel et al. 




10"' 



■ a " ■ ■ ■^■■■■i 



defined in lDaddi e t al. (2007b). 
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Figure 4. Plot of SFRt t/SFR{jy corr versus Lir for sBzK galaxies used 
for IR excess analysis. Filled squares are galaxies with spectroscopic red- 
shifts while open squares are galaxies with photometric redshift estimates 
only. The horizontal dashed line denotes the IR excess criterion; galaxies 
above the line are IR excess, below the line are IR non-excess. Black crosses 
represent galaxies that have been directly detected in the 4Ms CDFS and 
2Ms CDFN data. 



The z — K colour acts a proxy for stellar mass, probing the op- 
tical S EP about the Balmer and 4000 A breaks for galaxies wit h 
z ~ 2 jReddv et al.ll2005l : IShaplev et al]|2005l:IWuvts et aljEoolh . 
while the B — z colour corrects for reddening effects (see Equation 
[6). Additionally the sBzK selection was empirically calibra ted to 
make it independent of reddening effects dPaddi et al.ll2004l) . This 
is an important feature because emission from regions of intense 
star-formation is expected to be heavily reddened. Only galaxies 
that satisfy the sBzK criterion are subjected to 24^tm excess analy- 
sis. The BzK plot (see Figure[TJ also identifies passively evolving 
galaxies at z~2 (pBzK) and normal galaxies at z < 1.4 (nBzK), but 
these are not included in the subsequent IR Excess analysis. While 
the sBzK criterion primarily selects z ~ 2 galaxies there remains 
a fraction of low redshift interlopers (see Figure O. Additionally, 
a large spread in redshift distribution would also adversely affect 
the X-ray stacking analysis. Therefore the redshift range for sBzKs 
subjected to IR Excess analysis is narrowed to 1.2 < z < 3.0. The 
results of BzK selection and the subsequent IR Excess analysis are 
outlined in Table 1 . 

For the pur pose of thi s analy sis we calculate SFR/h following 
the relationship Kennicutt (1998): 



SFRmiM^yr' 1 ) = 1.73 x 10~ 10 L/h(L g 



(2) 



where Lir is the bolometric IR luminosity of the galaxy between 8 
and 1000 /J.m (Z/8-iooonm)- The observed 24/^m flux is k-corrected 
to rest-frame 8/xm flux u sing the IR star-form i ng gal axy templates 
of lCharv & Elbazl i feOOll) . The lCharv & Elbazl bOOlh templates are 
em pirically derived from o bservations of l ocal star-forming galax- 
ies. iMagnelli et all d2009h found that the ICharv & Elbazl 1 2001) 
templates exhibit a stronger correlation to the emission of ob- 
se rved galaxies at z ~ 1 th an the star-forming galaxy tem plates 
of iPale & Heioi] J2002l) and lLagache. Pole, & Pugetl J2003I) . The 
bolometric Lir is inferred from monochromatic 8^tm luminosity 
(i8nm) using the relationships for lCharv & Elbazl j200ll) templates 



4.31 



(3) 



if log( 



> 9.75 and 



log(- 



0.93log(^^-) + 1.23 (4) 



if log(^) <9.75. 

SFRc/y is estim ated using Equation [5] a relationship derived 
bv |Paddietal1 ( l2004l) : 

where L1500 is the 1500 A monochromatic luminosity. The UV 
spectrum is approximated to a power law derived from observed 
B,v,i and z bands (rest-frame UV at z ~ 2) and L1500 is then 
calculated. The dust reddening factor E(B — V) is estimated using 
the following criterion. 



E(B -V) = 0.25(B - z + 0.1)as 



(6) 



The UV emission is corrected for reddening, following the empiri - 
cally derived reddening correction method of lCalzetti et al.l (2000), 
and we use this reddening corrected UV emission to again calculate 
the SFR (SFRuv,corr)- To be defined as IR Excess sBzKs (here- 
after IRXs; Compton thick AGN candidates) a galaxy must satisfy 
Equation|7] 



log ( 



SFRw + SFRir , 



> 0.5 



(7) 



SFFLlJVcorr 

The sBzKs that do not satisfy this criterion are classified as IR Non- 
Excess sBzKs (hereafter IRNXs). 

If an apparent IRX (see Equation [Vj possesses a 24^im flux 
that is less than 3 times the background noise it is not considered 
a secure-detection and removed from the sample, as the measure- 
ment of an IR excess is not necessarily robust. Clearly it is neither 
necessary nor desirable to apply this cut to IRNXs. Between CPFS 
and CPFN there are 356 IRXs and 879 IRNXs. A K-S test demon- 
strates that the IRXs and IRNXs have significantly different redshift 
distributions, with null hypothesis p<0.001. The IRXs have higher 
redshifts on average (median z = 2.06) compared to the IRNXs 
(median z = 1.78), see Figure[3] The high density of spectroscopic 
observations at z — 1.6 is due to a cluster in the CPFS sample 
that has been extensively observed dCastellano et alj|2007r) . There 
is also a strong positive correlation between Lir and the ratio of 
SFRtot/ SFRuv,corr (the IR excess parameter) in Figure [4] As 
a consequence, the overwhelming majority of IR luminous sBzKs 
are classified as IRXs (85.3%; 304/356 with L m > 10 12 L Q ) as 
opposed to IRNXs (19.1%; 168/879 with Lir > 10 12 L Q ). 



2.3 X-ray data reduction 

The 4Ms CPFS Chandra data was obtained over the course of 3 
separate Chandra observing cycles: 1 Megasecond (Ms) between 
May 27th 2000 and Pecember 23rd 2000, another lMs between 
September 20th 2007 and November 4th 2007 and the final 2Ms 
between May 18th 2010 and July 22nd 2010. The 2Ms CPFN 
data was collected between 23rd February 2000 and February 22nd 
2002. All observations were taken using t he Advanced CCP Imag- 
ing Spectrometer imaging array (ACIS-I; lGarmire et al.ll2003l) . 
We now present a brief description of the data reduction 
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method which incorporates the use of CIAO data analysis soft- 
ware version 4.2 tools and c ustom software. Fo r a more detailed 
descrip tion of the method see lLaird et alj d2009h and lNandra et al.l 
din pred) . The raw data for each observation (henceforth, obsID) 
is corrected for known aspect offsets and hot pixels are removed. 
Then charge transfer inefficiency (CTI) and time dependant gain 
corrections are applied. Potential flare events are identified and re- 
moved in each obsI D using custom software (for a more details see 
Nandra et al. 20o3). We manually removed flares from 5 obsID's 
(441, 581, 1431, 2405 and 2406) in CDFS and 2 obsID's (966 and 
3389) in CDFN. 

An additional astrometric correction to each obsID is calcu- 
lated by comparing X-ray source positions to those of a refer- 
ence source catalogue with superior astrometry. Running the CIAO 
wavdetect tool on each obsID, bright X-ray sources are selected us- 
ing a threshold probability 10 -6 . These should be common to the 
majority of individual obsIDs and hence ideal for marrying up the 
astrometry. In CDFS we take reference astrometry from a catalogue 
of H-band selected galaxies (sources with the H-band magnitudes 
between 18 and 23 and no photometr ic flags) from the GOODS- 
MUSYC survey jGawiseretafl r2006). We removed the brightest 
sources to eliminate extended objects with potentially incorrect po- 
sitions, and faint magnitudes were omitted to reduce the likelihood 
of chance associations. The reference astrometry for 2Ms CDFN is 
from an r-band selected catalogue of ICapak et alj 12004J). We did 
not apply a magnitude cut to the r-band catalogue because it made 
no practical difference to the astrometry correction. 

Following the astrometry correction, merged images and ex- 
posure maps are created in four energy bands; 0.5-2 keV (soft 
band), 0.5-7 keV (full band), 2-7 keV (hard band) and 4-7 keV 
(ultra-hard band). Photons with energies < 0.5 keV and > 7 keV 
are ignored because the respons e of the ACIS-I instru ment drops 
off marked ly at these energie s {Weisskopf et al. 2002). Using the 
methods of lLaird et alj J2009I) , point source catalogues containing 
flux estimates and PSF's are generated for the 4Ms CDFS and 2Ms 
CDFN merged images. The full, hard and ultra-hard band fluxes are 
corrected to 0.5-10 keV, 2-10 keV and 5-10 keV respectively using 
a power-law model with Y — 1.4 and galactic Njj = 8.95 x 10 19 
cm -2 and N H = 1.51 x 10 20 cm -2 in CDFS and CDFN respec- 
tively. 

Our point source lists for 4Ms CDF S and 2Ms CDFN i mages 
were generated following the method of lLaird et alj d2009l) . find- 
ing 569 and 504 X-ray point sources respectively. The threshold 
poisson probability for sources in this catalogue is 4 x 10~ 6 (a 
more in depth de scription of the detection method can be found in 
lLaird et alj 2009) to ensure a low false identification rate. The 4Ms 
CDFS source catalogue o f IXue et alj J20llh contains 740 sources, 
of which 500 are coincident with sources from our catalogue when 
matched using a fixed 1.5 arcsecon d circular radius ( corresponding 
false match probability 1.0%). The lXue et alj feOllh catalogue has 
more point sources because of the lower threshold in their analysis. 
The 69 sources unique to our catalogue are low significance de- 
tections acco rding to ou r methodology. Likewise, the 240 sources 
unique to the lXue et alj J201 ll) catalogue are low significance de- 
tections according to their method, the majority of which have pois- 
son probabilities just below our detection threshold. 

2.4 Stacking Procedure 

We infer the average X-ray emission of the X-ray undetected sBzKs 
by stacking them. Stacking involves taking a sample of X-ray faint 
galaxies and summing their emission to create a long exposure im- 



age with an improved Signal-to-Noise (S/N) ratio. Data collected 
using Chandra ACIS-I can be summed to creat e extremely long ex - 
posures without reaching the confusion limit dBrandt et alJl200lT) . 
so is ideal for stacking analys is. The stacking tech niqu e employed 
here i s based on the work of iNandra et alj (2002) and lLaird et alj 
(2005). Fixed diameter circular apertures are used (as opposed to 
variable aperture techniques) to extract the source counts, as this 
has proved the most effecti ve and reliable so urce count extraction 
method for X-ray stacking (Laird et al. 2005). A fixed aperture of 
radius 1 .5 arcseconds provides the best S/N ratio in both 4Ms CDFS 
and 2Ms CDFN. Sources at large off axis angles (OAA) have wider 
PSFs, therefore the source counts become more diffuse and less 
significant relative to the background. Thus, including sources with 
large OAA's can have a negative effect on the S/N ratio of the final 
stacked image. To overcome this problem we define a maximum 
OAA, known as the 'inclusion radius', for sources that are to be 
stacked. We find an inclusion radius of 5.5 arcminutes gives the 
optimum S/N ratio for stacks in both CDFS and CDFN. 



Prior to stacking we mask out any sources that are within 
an area defined by the 2 times 90% PSF of an X-ray detected 
source. This prevents bright sources from contaminating the stack 
and skewing the low count statistics of the faint sources we are in- 
terested in. A local background is estimated for each source by ex- 
tracting counts in an annulus centred on the candidate source with 
inner and outer radii of 10 arcseconds and 30 arcseconds respec- 
tively. Pixels from which the background counts are extracted are 
chosen by randomly shuffling the position 10,000 times within the 
annulus and averaging the counts per pixel. Any part of the annu- 
lus that overlaps with the 2 times 90% PSF of an X-ray detected 
source is masked out to prevent bright source contamination. X- 
ray undetected candidate sources within the background annulus 
are not masked out because they are not confirmed X-ray sources 
and must still be treated as part of the background. The total back- 
ground counts are scaled to match the source extraction aperture in 
both area and exposure. 
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Table 2. X-ray properties of sBzK galaxies detected in the 4Ms CDFS. Column (1): X-ray source name, galaxies suffixed with f have been detected in the 
lMs CDFS also; column (2): BzK subsample; columns (3) and (4): right ascension (RA) and declination (Dec.) of source; column (5): redshift, redshift's are 
spectroscopic where available (suffixed with s ), otherwise photometric redshifts are used (suffixed with p ); column (6): soft band flux, units 10~ 16 erg cm -2 
s _1 ; column (7): hard band flux, units 10 -16 erg cm~ 2 s _1 ; column (8): rest-frame hard band Luminosity derived from the soft band flux, units 10 42 ergs 
s _1 ; column (9): HR, HR= (H — S)/(H + S), where H and S are hard- and soft-band count rates respectively. 



Name 


Tvrje 


RA 


Dec 


Redshift 


1 U.o — Zke V 


1 2. — 1 U/c e V 


^ Z— LUkc V 


HR 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


cdfs4Ms.031t 


IRX 


53.15731 


-27.87007 


1.60 s 


65.69tii| 


119.92±i;« 


88.47 + i| 9 


— n 4Q+0.01 

U - 4C> -0.01 


cdfs4Ms.043t 


IRX 


53.16146 


-27.85597 


2.91P 


0.841°;^ 


7 ,1+1.09 
' - ol -0.90 


7.78l 1 1 ; 4 8 


qn+0-l° 


cdfs4Ms.077t 


IRX 


53.14879 


-27.82116 


2.58 s 


0.491°;^ 


=9+0.69 
°- oz -0.62 


9 Q9+0.65 


91 +013 
u - zl -0.12 


cdfs4Ms.078t 


IRX 


53.1801 


-27.82061 


1.92 s 


28.nt°,; 7 T 5 2 


58.78l|j? 


66.15li;^ 


n o 7 +0.02 
u - a '-0.02 


cdfs4Ms.079 


IRX 


53.05903 


-27.81944 


2.28P 


n 2 =+0.09 


< 1.10 


1 02+ - 36 
i - UzI -0.31 


-0.08l°; 3 22 


cdfs4Ms.091t 


IRX 


53.18579 


-27.80996 


2.59 s 


3.76l°; 29 


11.1211;" 


22.68li;S 


-0-2il°;° D 


cdfs4Ms_120t 


IRX 


53.18344 


-27.77658 


2.69 s 


5.56t°:SS 


8-04l i° 


37.53l 2 ; 3 [ 


-0 5i+o- 05 

u - oi -0.05 


cdfs4Ms_123t 


IRX 


53.04907 


-27.7745 


2.21P 


Q = + 0.45 

°- oz -0.40 


25 24 +1 - 74 


31.0911:2 


-0 91+0.04 
u - zi -0.04 


cdfs4Ms_126t 


IRX 


53.13114 


-27.77308 


2.22 s 


o-6o±g:i? 


i-92±S:SS 


9 99+0.44 

z - zz -0.40 


-n 14+0-1 6 


cdfs4Ms_131t 


IRX 


53.09394 


-27.76773 


1.73 s 




35.16li;H 


c = =+0.46 
O-OO-0.39 


n 4o+ - 04 

u - 4u -0.04 


cdfs4Ms_144t 


IRX 


53.12491 


-27.75832 


1.22 s 


52.01_ 99 


124.08!!; 44 


30.26+81? 


-0 S1+ - 01 


cdfs4Ms_161t 


IRX 


53.04549 


-27.73748 


1.62 s 


54.79±i°| 


114.69ti;g 


75.51±i : S 


-0 S7+°' 01 
u - a '-0.02 


cdfs4Ms_175t 


IRX 


53.10697 


-27.71826 


2.29 s 


13.07±g;|f 


55.59 + 2-|g 


53.40 2 j 26 09 


-n 04+ ' 03 

u - u4 -0.03 


cdfs4Ms_183t 


IRX 


53.1337 


-27.69866 


1.39 s 


6-68l ; 039 


i6.35i;; 9 7 ° 2 


= ok + 0.38 
0-00-0.34 


-0 ,n +0 '° 5 


cdfs4Ms.215t 


IRX 


53.15812 


-27.88552 


2.73P 


o-9i±S:ii 


< 1.34 


6.44li' 2 ? 


-0 49+ - 24 

u -* a -0.24 


cdfs4Ms.224 


IRX 


53.17406 


-27.85968 


1.93 p 


0.3818:15 


l.08l°; 6 o 2 8 


Q1+ ' 31 
u - al -0.28 


n nc + 0.29 
-0.06Io.22 


cdfs4Ms.264f 


IRX 


53.05881 


-27.70842 


2.03 s 


1 62+ ' 28 


7-57±l;™ 


4-5118:66 


02+°- 13 
u - uz -0.12 


cdfs4Ms.334 


IRX 


53.16222 


-27.71213 


2.23P 


u ' ou -0.19 


6.56lJ-|| 


l.89l°; 7 7 6 3 


n 42+o-i 5 

U - 4Z -0.15 


cdfs4Ms_369 


IRX 


53.06313 


-27.69953 


2.40 s 


l.76l°; 3 2 ° 5 


4.2511;^ 


°* oz -1.20 


-0.2518;!! 


cdfs4Ms.434 


IRX 


53.19878 


-27.84391 


1.51 p 


0.281°;" 


8.6411:3* 


0-3118: 1 ! 4 


0-69l°;° 9 


cdfs4Ms.438 


IRX 


53.16685 


-27.79876 


2.00 s 


26+ 09 
u - ZD -0.08 


0.94l°|0 


0.68l°; 24 


-o.02l°; 2 I 


cdfs4Ms.451 


IRX 


53.1312 


-27.84129 


1.61 s 


< 0.08 


i-06l8:5? 


< 0.11 


0.48l°; 2 l 


cdfs4Ms.460 


IRX 


53.09755 


-27.71552 


2.14 s 


< 0.14 


4.4011; 29 


< 0.46 


0-76l°; 2 9 


cdfs4Ms.515 


IRX 


53.05683 


-27.79852 


1.83P 


n 14+0.08 

-0.07 


< 0.60 


0.28l°;l 6 5 


—0 1Q+ ' 42 
u - la -0.35 


cdfs4Ms.013t 


IRNX 


53.19605 


-27.89266 


2.74P 


= 4 = +0.40 


22.19l 2 -°° 


39.08 + 2-g| 


n nc + 0.06 
-°- 06 -0.05 


cdfs4Ms.042t 


IRNX 


53.15067 


-27.85735 


1.61 s 


9 i =+0.23 
z - io -0.19 


fi 12 +0.97 
D - ±z -0.80 


2.95t8-iJ 


n 99 + O.O8 
u - zz -0.08 


cdfs4Ms.050t 


IRNX 


53.07463 


-27.84865 


1.54 s 


1 co+0.19 
l-t>»_0.18 


47.66l 2 ; 2 3 


1.87°_- 23 21 


n 7Q+0-03 


cdfs4Ms.083t 


IRNX 


53.04569 


-27.81555 


1.39 p 


11 AO+0. 48 
ii.ty_ 4g 


22.20ll;| 4 


9-99l°; 42 


-o 4n+°'° 3 

u - 4u -0.03 


cdfs4Ms_088t 


IRNX 


53.14987 


-27.81405 


1.31 s 


i 99 + O.I6 

-■-■^^ — 0.15 


, go+0.56 
1 - MO -0.52 


0.8818;! 2 


-n 44+ - 11 

u -*^-0.11 


cdfs4Ms.089t 


IRNX 


53.17929 


-27.81254 


1.73 s 


1 99+O.I7 
J - 4 - 4 — 0.16 


19.82±l;|f 


9 n7 +0.29 
z - u '-0.26 


n ==+0.06 

u.dd_ 05 


cdfs4Ms_133t 


IRNX 


53.16283 


-27.76716 


1.22 s 


16.02l°| 9 


52.66 + 2-?4 


o 27 +0.34 
y - z ' -0.32 


n 90+0.03 

u - zo -0.03 


cdfs4Ms_155t 


IRNX 


53.02409 


-27.74643 


1.61 s 


o 7O+0.33 
J -< 8 -0.29 


S.Olll; 49 


5-141811 


-0 S5+°'° 9 
u - oo -0.07 


cdfs4Ms_157t 


IRNX 


53.16272 


-27.74424 


2.52P 


= ,1+0.35 
°" L1 -0.31 


20.131^3 


28.1311;? 2 


-0 O7 +0 '° 5 
u - u ' -0.05 


cdfs4Ms.206t 


IRNX 


53.08925 


-27.93046 


2.61 s 


c. n 7+0.42 
°- u ' -0.38 


14.161 2 '? 2 


31.25l 2 ; 6 33 


n 1 0+0. 07 

-U. 18-0.08 


cdfs4Ms.208t 


IRNX 


53.20484 


-27.91799 


2.02P 


7 91 +0.52 
' - z± -0.47 


13.50±l;f| 


19.87li; 4 3 4 


_ n 41+0.07 

u - 4i -0.08 


cdfs4Ms.226t 


IRNX 


53.0601 


-27.8529 


1.54 s 


n 4 =+o.i3 

u '^°-0.12 


6-83l ; 94 


n ck + 0.16 
u - 00 -0.14 


50+ - 12 
u - ou -0.09 


cdfs4Ms.246 


IRNX 


53.18171 


-27.78293 


1.57 s 


o-3il°o; 10 


< 0.71 


QO+0.14 


-0.37 9 j 33 32 


cdfs4Ms.258 


IRNX 


53.06154 


-27.73409 


1.67 p 


1.081°;" 


O Ol +1.16 


1 65+ - 29 
l.OO_ 27 


-0-i6l8:i o 


rHf«4M« 971 


IRNX 




77 701 1 Q 


1 9^ s 


n 55+ 023 

U.OO_g 22 


9 97+I.9O 
-1.03 


n oo+0.13 
U - d,J -0.13 


O5+0- 29 
U.UO_ 22 


cdfs4Ms_283 


IRNX 


53.08729 


-27.92957 


2.88^ 


1 qq+0-32 

l.»0_Q 2g 


18.48l 2 -l 9 


16.191 2 ] 2 


o = +0.08 


cdfs4Ms.321 


IRNX 


53.02796 


-27.74872 


2.38P 


n 70+O.I8 
u - '^-0.17 


1 QQ+i-iO 
1-99-0.94 


07+0. 85 
•J- 31 -0.79 


-0 12+ - 26 

u - lz -0.19 


cdfs4Ms.370 


IRNX 


53.05303 


-27.69705 


1.78 p 


1 ci+0.31 
± - o± -0.26 


4.06+i' 92 


9 01 +0.57 
2-81-0.49 


-0.1918;?! 


cdfs4Ms.405 


IRNX 


53.041 


-27.83608 


1.78 p 


20+ ' 11 

u - zu -o.io 


1 94+ - 84 


n 00+O.2O 


0.28l8; 2 « 


cdfs4Ms_441 


IRNX 


53.09494 


-27.75793 


1.87 p 


< 0.07 


1 fif) +0.66 
1- DU _0.62 


< 0.15 


0.67l°; 2 ^ 


cdfs4Ms_459 


IRNX 


53.12284 


-27.72285 


1.63 p 


< 0.14 


3 85+ 1 ' 18 


< 0.20 


70+0- 23 
u - ' u -0.12 


cdfs4Ms.494 


IRNX 


53.1255 


-27.88646 


2.64 p 


o.66±°:iS 


< 1.11 


A 1O+0.95 
4.18_ 087 


_n 46+ - 26 

u - 40 -0.26 


cdfs4Ms.498 


IRNX 


53.12022 


-27.79886 


1.38 s 


n 90+0.10 

U.2»_ 08 


< 0.48 


o.24±8:8? 


-0 5n+°' 26 

U.OU_ 32 


cdfs4Ms.503 


IRNX 


53.16382 


-27.7735 


2.70 p 


15+ - 09 


< 0.66 


i-05l8: 5 5 9 3 


-0 17+°- 4 l 
u . 1( -0.33 


cdfs4Ms.512 


IRNX 


53.11056 


-27.8236 


1.47 s 


24+°' 09 
u - z4 -0.07 


< 0.32 


25+°-° 9 
u - ZJ -0.07 


-0 7n+ 08 

u -' u -0.30 


cdfs4Ms.301 


pBzK 


53.24571 


-27.86106 


1.10 p 


3.02l°|? 


9-08li-.II 


1 99+O.I4 
1-^^-0.13 


n iR+O-ii 
-0.18_ 010 


cdfs4Ms.318 


pBzK 


53.04476 


-27.77442 


1.61 s 


o.8o±°:iS 


< 1.29 


i-08l°; 2 9 


_n 45+ - 22 

u - 4 °-0.20 


cdfs4Ms.457 


pBzK 


53.02676 


-27.76529 


1.33 s 


< 0.21 


ok + 117 
z - yd -1.09 


< 0.16 


4Q+0- 25 
u - 4y -0.19 
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3 RESULTS 

3.1 Direct Detections 

We match the sBzK galaxies in 4Ms CDFS and 2Ms CDFN to X- 
ray point sources using a 1.25 arcsecond search radius to establish 
direct detections. The BzK source density is sufficiently low that 
at 1.25 arcseconds all matches are unique, and the astrometry cor- 
rections detailed in Section 1231 have made the overall offsets be- 
tween the optical and X-ray positions relatively small. We directly 
detect 5.8±0.8% (49/851) of sBzKs in the 4Ms CDFS image and 
5.4±1.0% (32/594) of sBzKs in the 2Ms CDFN image (see Tables 
2 and 3). The rest-frame 2-10 keV luminosities (Z/2-iofcev; derived 
from F . 5 -2kev) of sBzKs range from 1.1 x 10 41 — 3.5 x 10 44 erg 
s" 1 with median Z/2-iofceV = 3.5 x 10 42 erg s _1 . Approximately 
70% (57/81) of the detected sBzKs possess L 2 -iokeV > 10 42 erg 
s" 1 and it is reasonable to assume all of these objects host AGN. 
The X-ray emission of the remaining 26 sBzKs is probably also 
mostly from AGN, but a significant contribution from stellar pro- 
cesses is also possible for these lower luminosity X-ray sources. 

The detection rate of X-ray sources is a strong function of 
Lir (see Figure . This agrees with previous observations that 
AGN are more prevalent amongst the most IR luminous galaxies 
dSanders & Mirabelll99fj|Risaliti et al.l200d) . The HR distribution 
of X-ray detected sBzK galaxies is notionally bimodal (Figuref5](, 
covering a wide range of HR (-0.70 to 0.82). Using a simple photo- 
electric absorption model, an AGN at z ~ 2 with power-law emis- 
sion r = 1.9 obscured by material with Nh > 10 23 cm -2 should 
possess a HR > -0.08. Therefore 44% (36/81) of the X-ray de- 
tected sBzKs possess Nh = 10 23 cm -2 or greater. This suggests 
the sBzKs are a mixture of obscured AGN, unobscured AGN and 
purely star-forming galaxies. 

Splitting the sBzK galaxies into IRX and IRNX sub-samples 
we directly detect 13.2±1.9% (47/356) and 3.9±0.7% (34/879) re- 
spectively. It is clear that a greater proportion of IRXs than IRNXs 
are detected in X-rays. The IRXs are also significantly more X-ray 
luminous (median Z/2-iofceV = 6.44 x 10 42 erg s" 1 ) than IRNXs 
(median L2— lofcev = 2.26 x 10 42 erg s" 1 ). On the other hand, 
the higher detection fraction of IR Excess sBzKs may be more re- 
lated to their IR luminosity, rather than the IR excess per se. Fig. 4 
clearly demonstrates a very strong correlation between the excess 
parameter and Lir, such that almost all of the most IR-luminous 
BzK galaxies also exhibit an IR excess. 

If we bin the sBzKs by IR luminosity and then compare IRXs 
and IRNXs, as in Figure[6] the detected fractions are largely similar. 
If anything, when controlled for IR luminosity, the IRNXs show a 
higher detected fraction, although this is not statistically significant. 
A similar fraction of IRXs and IRNXs possess Z/2-iofceV > 10 42 
erg s _1 (74% and 65% respectively), so the conclusion should 
also hold for the AGN fraction as well (i.e. as opposed to the 
detected fraction), especially when considering that the strongest 
star-formation will be in the most IR-luminous (and hence excess) 
galaxies. 

The IRXs and IRNXs have very similar distributions of HR 
(see Figure |5). There are 14 (30%) IRXs and 1 1 (32%) IRNXs 
with HR > 0.2. These hard X-ray sources are almost certainly 
obscured AGN candidates, but there is no evidence that there is 
a heavier concentration of heavily obscured AGN within the IRX 
sample compared to the IRNX sample. A K-S test shows that, 
based on their HR distributions, the IRXs and IRNXs are consis- 
tent with each other (null hypothesis probability p=0.829). There 
are 4 sBzKs in our sample in CDFS that have previously been 
confirmed using X-ray spectroscopy as containing Compton thick 
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Figure 5. Histogram of HR of directly detected sBzK galaxies. The black 
line represents all sBzKs, the red line represents IRXs and the blue line 
represents IRNXs. 
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Figure 6. Plot of detected fraction of sBzKs against Lju where red di- 
amonds are IRXs, blue diamonds are IRNXs and black diamonds are all 
sBzKs. Mult inomial errors for s BzK detected fraction calculated using re- 
lationship of Nandra et al. I 20q3). 



AGN: cdfs4Ms_077 and cdfs4Ms_28 3 by iFeruglio et aljfcoilh 
and cdfs4Ms_126 and cdfs4Ms_226 bv lBrightman & Uedj d2012bl) 
(note t hat cdf s4Ms_077 was also identified by iBrightman & Ueda 
l2012bl as containing a Compton thick AGN). Only cdfs4Ms_226 
is an IRNX, while cdfs4Ms_077, cdfs4Ms_126 and cdfs4Ms .283 
are IRXs. The fraction of IRXs confirmed as Compton thick 
AGN from their X-ray spectra is lower than expected (3/49), 
but the X-ray detected IRX sample is too small to draw any 
strong conclusions regarding the nature of IRXs as a whole. 



3.2 Stacking sBzKs 

In addition to analyzing the direct detections we also stacked the X- 
ray undetected sBzKs in 4Ms CDFS and 2Ms CDFN to probe their 
average X-ray properties. In the first instance the X-ray photons 
are stacked in 3 broad bands: 0.5 — 2 keV (soft), 2 — 7 keV (hard) 
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Table 3. X-ray properties of sBzK galaxies detected in the 2Ms CDFN. Column (1): X-ray source name; column (2): BzK subsample; columns (3) and 
(4): right ascension (RA) and declination (Dec.) of source; column (5): column (5): redshift, redshift's are spectroscopic where available (suffixed with s ), 
otherwise photometric redshifts are used (suffixed with p ); column (6): soft band flux, units 10 -16 erg cm -2 s _1 ; column (7): hard band flux, units 10 — 16 
erg cm -2 s — 4 ; column (8): rest-frame hard band Luminosity derived from the soft band flux, units 10 42 ergs s — 1 ; column (9): HR, HR= (H — S) /(H + S), 
where H and S are hard- and soft-band count rates respectively. 



NtllTIC 


T 

type 


RA 


Dec 


Redshift 


-^0.5 — 2keV 


-^2 — lOfceV 


r 

^2 — 10fceV 


HR 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


hdfn.015 


IRX 


189.13587 


62.13338 


1.99 s 


1 qrj + 0.31 

J - OZ -0.28 


13.991?;!! 


o 40+0.81 
•J-^y-0.74 


0.3518:1? 


hdfn.050 


IRX 


189.06689 


62.18556 


1.91 p 


1 V K9 + 0.85 

i(.oz_ 77 


149.24±|;!§ 


40.57li'^ 


o-28l° ;™ 


hdfn.060 


IRX 


189.22245 


62.19449 


1.27 s 


n 77+0.19 

u -' '-0.16 


< 1.25 


f) n +0.12 
U - Oi -0.11 


— n 47+0- 20 

u - 4 ' -0.20 


hdfn.067 


IRX 


189.23252 


62.20028 


2.74 s 


1 Qn+0' 23 


1 r;o + 0.80 

x.oo_ 73 


q 09+1. 64 


— n 59+0.14 

u - oz -0.15 


hdfn.080 


IRX 


189.14395 


62.2115 


1.22 s 


n 7C+O.I8 


z - 44 -0.74 


0.45 +0 o 9 n 


-n 16+° 18 

U - 1D -0.17 


hdfn.082 


IRX 


189.26096 


62.21223 


2.31 p 


9 4O+0.31 
z - 4y -0.28 




10.46li;f 9 


— n 45+ ' 09 
u - 40 -0.11 


hdfn.098 


IRX 


189.35108 


62.23327 


1.59 p 


K OQ+0' 51 
J - za -0.43 


55.17ti;« 


6.95ir 5 6 


n ^7+0.04 


hdfn.107 


IRX 


189.14843 


62.24015 


2.01 s 


1 7O+0-26 
J- '3-0.24 


20.20±f°l 


4.8ilg:S 


n Qo+o.07 


hdfn.l 19 


IRX 


189.21598 


62.25136 


2.18 p 


7 in+0.55 
'• lu _0 .48 


22.20±i: 2 6 


24.82+1;?? 


-n 21 + - 05 

u - zl -0.06 


hdfn.120 


IRX 


189.25198 


62.25248 


2.29 p 


01 +0.38 
°- ol -0.36 






-n n7+ ' 08 

u - u '-0.09 


hdfn.127 


IRX 


189.34672 


62.2607 


2.24 s 


4 42 +0.47 


15.77t 2 ;« 


16.90±i;gl 


-n i4+o- 09 

u ' 14 -0.07 


hdfn.137 


IRX 


189.27041 


62.26714 


1.88 p 


27.65l ;° 4 T 


80.18±|;g| 


6O.95l2.14 


— n 95+0. 03 

u - zo -0.03 


hdfn.144 


IRX 


189.31159 


62.27141 


1.52 s 


s.93l8:S 


24.56t2- 4 g 


in 24+ ' 72 


n 97 +0.06 
—U- £ < -0.05 


hdfn.158 


IRX 


189.17957 


62.28656 


2.28P 


"■i7lo:S 


30.20+2 


44 Q4+ 2 - 94 


-n.28l° : ° 


hdfn.168 


IRX 


189.42735 


62.30341 
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and 0.5 — 7 keV (full). The stacked photon count rates are then 
converted to fluxes assuming a V = 1.4 power-law: Fo.s-2keV 
(soft), F 2 -iokeV (hard), and F .5_i fceV (full). Starting with the 
full sample, we stacked 370 sBzKs in 4Ms CDFS and find a highly 
significant signal in the soft and hard band. The stacked signal pos- 
sessed a HR = -0.53 ± 0.08. We stacked 318 sBzKs in 2Ms 
CDFN and again found significant soft and hard band signals, with 
a HR = -0.47 ± 0.09. The stacked flux of the X-ray undetected 
sBzKs in CDFN is higher than for those in CDFS (see Table 5) be- 
cause the 4Ms CDFS is deeper, and the brightest sources stacked 
in the CDFN 2Ms image would be detected and resolved out in 
the CDFS 4Ms image. For the total sample, both CDFS and CDFN 
sBzK galaxies exhibit soft HRs in these broad bands and there is no 
evidence to suggest a substantial proportion of the sBzKs contain 
heavily obscured AGN. 

Splitting the sBzKs using the IR excess criterion (Equation 
0, we stacked X-ray undetected 72 IRXs and 298 IRNXs in 4Ms 
CDFS. The stacked signal of the X-ray undetected IRXs is signifi- 
cantly detected in both soft and hard band with a HR = —0.48 ± 



0.11. The stacked signal of X-ray undetected IRNXs is also signif- 
icantly detected in the soft and hard band, but is more X-ray faint 
(see Table 5) and has a softer HR (—0.55 ± 0.11). The higher flux 
of the stacked IRXs is in agreement with our observations of di- 
rectly detected sources (Section [3.1t . Again, however, this can be 
attributed to the higher IR luminosity for the excess sources, whose 
median Lir = 3.6 x 10 12 ' Lq is over an order of magnitude higher 
than the IRNX sources (median Lib = 3.1 x 10 11 L Q ). Addition- 
ally the HR and effective T of X-ray undetected IRXs and IRNXs 
in CDFS are consistent with each other at the la level (see Table 5) 
which suggests the composition of the sub-samples is similar. 

Turning to the 2Ms CDFN we stacked 107 X-ray undetected 
IRXs and 211 X-ray undetected IRNXs. As with the stacking in 
4Ms CDFS, the IRXs are significantly detected in both the soft and 
hard bands. The HR (— 0.50±0.12), agrees with the stacked signal 
of 4Ms CDFS X-ray undetected IRXs at the la level. The X-ray 
undetected IRNXs are significantly detected in the soft band (S/N 
= 8.2) but the hard band is a low significance detection (S/N = 2.7). 
The stacked fluxes are again higher than 4Ms CDFS X-ray unde- 
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tected IRNXs (see Table 5) and the HR = -0.45 ± 0.15. The HRs 
of X-ray undetected IRXs and IRNXs are consistent at the la level 
although the low significance hard band detection of X-ray unde- 
tected IRNXs means the HR is probably softer than our current esti- 
mate. Ultimately all stacks of X-ray undetected sBzKs and the IRX 
and IRNX subsamples agree at the la level (Table [5} and there is 
no evidence from broad band stacking that heavily obscured AGN 
are more prevalent in IRXs than IRNXs. 

As discussed in Section |2~2l the IR excess criterion (Equation 
is positively correlated with redshift and Lir (Figures [3] and [4] 
respectively) and this is likely to account for the higher X-ray fluxes 
of the IRXs. To investigate this further, we have stacked samples 
binned by redshift and Lir. The sBzK galaxy samples from CDFS 
and CDFN have been combined to maximise the S/N ratio in each 
redshift and Lir bin. 

The X-ray undetected sBzKs are split into 3 redshift bins with 
ranges that have been chosen to distribute the sources as evenly as 
possible: 1.2 < z < 1.8, 1.8 < z < 2.4 and 2.4 < z < 3.0. 
The stacked emission of X-ray undetected sBzKs with 1.2 < z < 
1.8 and 1.8 < z < 2.4 possess HR = -0.47 ± 0.10 and HR= 
—0.48 ± 0.09 respectively and exhibit similar X-ray fluxes (see 
Table 6). The effective T values are consistent with those observed 
for X-ray undetected IRX and IRNX subsamples at the la level. 
The 2.4 < z < 3.0 X-ray undetected sBzKs have a similar flux 
to the lower redshift sources in the soft band, but no detection in 
the hard band (S/N = 0.6). The HR is therefore very uncertain, but 
inferring the hard flux from the full band detection shows X-ray 
colours which are entirely consistent with the lower redshift bins. 
As the fluxes of each redshift-binned subsample are all similar, this 
naturally implies that the higher redshift sources are on average 
more luminous in the X-ray, but this is likely to be a selection effect 
due to the fact that the parent sample is magnitude limited, meaning 
that the higher redshift sources need to be more luminous overall 
to enter the sample. 

The X-ray undetected sBzKs were also split into 4 Lir bins 
to explore the connection between Lir and the IR excess crite- 
rion (Section^: L IR < 10 U L©, IO^Lq < Lir < 1O 12 L , 
1O 12 L < Lir < 10 13 L Q and Lir > 1O 13 L . There is a strong 
correlation between X-ray brightness and Lir of the sBzKs (see 
Table 6). This result is in agreement with our observations of di- 
rectly detected X-ray sources. All bins are significantly detected in 
hard and soft bands except for Lir < 10 11 Lq which has a low 
significance detection in the hard band (S/N=2.1). The HR values 
of all the Lir binned sBzK sub-samples are consistent with one 
another at the la level and the stacked signal of IRXs and IRNXs 
in CDFN and CDFS. As with the redshift selected galaxies, none 
of the Lir bins appear to preferentially select heavily obscured 
AGN, but the BzK galaxies which are more luminous in the IR are 
brighter in the X-ray. 



3.3 Stacking Smaller Bands 

When stacked in broad bands (Section [3.2b the X-ray undetected 
IRXs and IRNXs both appear to have relatively soft X-ray spec- 
tra, with effective T ~ 2. Based on this analysis there is no ev- 
idence for obscured AGN in either subsample. The broad bands, 
however, may mask more complex spectral features and give a mis- 
leading picture of th e actual spectral shape and h ence the level of 
obscuration (see e.g. Brightman & Nandra 2012al) . We have there- 
fore stacked the sBzKs in narrower energy bands to obtain a more 
accurate representation of their spectra (Table 4 and Figure|7J. The 
bands in Table 4 have been designed to provide the optimum bal- 
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Figure 7. Flux estimates for narrow band stacking of IRXs (red diamonds) 
and IRNXs (blue diamonds) galaxies in the 4Ms CDFS and 2Ms CDFN. 
The horizontal error bars represent the width of the energy band that is 
being stacked. The IRXs stacked in the two highest energy bands are not 
significantly detected. 



ance of resolution and signal strength. Narrower bands have fewer 
counts so the CDFS and CDFN sBzKs have been combined to max- 
imise the S/N ratio in these bands. This analysis shows clearly that 
the naive interpretation of the broad band colours gives a mislead- 
ing impression of the true spectral shape, which does not resemble 
a single power law. The total spectrum is clearly "composite" in 
nature with a soft component dominating the soft energies, and a 
flat, hard component coming through at high energies. The latter is 
a clear signature of obscured AGN in the sBzK galaxies. The soft 
component is more ambiguous and could originate from a range of 
sources. There could be contributions from unobscured AGN, scat- 
tered light in heavily obscured AGN and/or X-rays from pure star 
forming processes. The stacked spectrum indicates that there is a 
mixture of these effects within the sample. 

Splitting the sample into IRXs and IRNXs, we observe that 
they have remarkably similar spectral shape with the only notice- 
able difference being that IRNXs are consistently less X-ray bright 
(see Figure [7), just as was observed for broad band stacking (see 
Section [3~2l . Thus, while the narrow band analysis confirms the 
presence of residual AGN activity in the population of sBzK galax- 
ies as a whole, there is no evidence that obscured AGN or AGN 
generally are preferentially hosted by those with excess emission 
at 24 /jm. 



3.4 Passive galaxies 

There are 21 passively evolving BzK's (pBzK) in CDFS and 5 
pBzK's in CDFN. We directly detect 3 pBzK's in 4Ms CDFS and 
1 pBzK in 2Ms CDFN, see Tables 2 and 3. While these num- 
bers are small, and hence the uncertainties large, the pBzK detec- 
tion fraction (15±8.0%; 4/26) is actually greater than for sBzKs 
(7±0.7%; 81/1235). The sources cdfs4Ms_301, cdfs4Ms_457 and 
hdfn_045 have hard signals which suggest the presence of ob- 
scured AGN. The source cdfs4Ms_318 lacks a hard band de- 
tection but L0.s-2.0keV = 1-1 x 10 42 erg s _1 (derived from 
Fo.5-2kev)- The X-ray emission probably originates from ei- 
ther unobscured AGN activity or from a hot gas halo within 
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Table 4. Narrow band stacking of X-ray undetected sBzK galaxies in 4Ms CDFS and 2Ms CDFN. Column (1): energy range of narrow band; column (2): 
number of sBzKs stacked; column (3): S/N of stacked sBzKs; column (4): narrow band flux of stacked sBzKs, units 10 — 18 erg cm™ 2 s" 1 ; column (5): 
number of IRXs stacked; column (6): S/N of stacked IRXs; column (7): narrow band flux of stacked IRXs, units 10~ 18 erg cm™ 2 s~ 1 ; column (8): number 
of IRNXs stacked; column (9): S/N of stacked IRNXs; column (10): narrow band flux of stacked IRNXs, units 10 -18 erg cm -2 s _1 . 
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the galaxy dForman. Jones. & TuckeJll985l : [Sarazinlll997h . Intense 
star-formation as the origin of the observed X-ray emission is ex- 
tremely unlikely, recalling that these objects have been selected as 
passively evolving galaxies. Stacking the X-ray emission of the 13 
undetected pBzKs in CDFS and CDFN yields a marginal signal in 
the soft band, with average fo.5-2fceV = 2.23 ± 0.62 x 10 _1 erg 
cm" 2 s - . Stacking the hard band flux does not result in a signifi- 
cant detection. 



4 COMPARISON TO PREVIOUS RESULTS 

Daddi et al. | d2007bl) stacked the X-ray emission of 59 IRXs and 
175 IRNXs that were not directly detected in X-rays using a lMs 
CDFS image and source list. Their sample of IRXs had a much 
harder stacked signal (effective T = O.Stg'*) than our stacked sam- 
ple in the 4Ms CDFS and 2Ms CDFN data (see Table©. There are 
a number of possible reasons for this discrepancy, including differ- 
ences in the parent BzK sample and deeper X-ray data in our case. 
To help understand this discrepancy we stacked our samples of 
IRXs and IRNXs using only the lMs CDFS images and source list 
to provide a direct comparison (Table[7}. Due to the shallower data, 
fewer sBzKs are directly detected in the lMs CDFS so we stack 
more X-ray undetected IRXs (80) and IRNXs (311). Our samples 
of X-ray und etected IRXs and IR NXs are different to (and larger 
than) those of lDaddi et al] (2007b) because of the updated photom- 
etry of GO ODS-MUSIC v2 in cluding deeper 24/^m data used in 
our work dde Santis et al.lr2007h . 

The X-ray undetected IRXs were significantly detected in 
both the soft (S/N = 5.6) and hard (S/N = 3.6) bands with HR 
= — 0.16±0.16 and corresponding effective F = I.SIq'j. This ef- 
fective r is clearly rather harder the n we obtained earlier, and while 
formally softer than that derived by Daddi et al. ( 2007^) it is con- 
sistent with their results of at the la level. The X-ray undetected 
IRNXs are significantly detected in the soft band (S/N = 7.6), but 
the ha rd band is a lo w significance detection (S/N = 2.2), something 
which lDaddietal.1 d2007bl) also found. The HR = -0.52 ± 0.17 
and corresponding effective T — 2.lt '^ are poorly constrained 
but are also consistent with the effective F ~ 1.7 of lDaddi et al.l 
d2007bl) at the la level. 

Our analysis, and that of lDaddi et al.l d2007bl) . shows that the 
spectrum is more complex than a simple power, so a naive com- 
parison of the two spectral indices does not give the full picture. 
Nonetheless it is clear that we find that the overall spectrum of X- 
ray undetected IRXs is softer. The above analysis shows that this 



can be partially, but not completely explained by the difference in 
the parent BzK sample. A further possibility is that the lMs stack 
is contaminated by hard sources undetected in the lMs image, but 
detected and resolved out in the 4Ms data. There are 22 sBzKs 
(9 IRXs; 13 IRNXs) which are detected in the 4Ms CDFS image 
but are undetected (i.e. fall below the detection limit) in the lMs 
CDFS image (see Table |23). Of these 22 sBzKs, 6 IRXs and 6 
IRNXs were stacked in the lMs CDFS image. The 3 IRXs and 6 
IRNXs that were directly detected in the 4Ms CDFS images were 
omitted from the stacking as they fell outside the 5.5 arcminute in- 
clusion radius (see Section [2~4"t . Focusing on the X-ray detected 
IRXs (because their signal softened more dramatically) the aver- 
age HR = 0.14 ± 0.09, with effective T = 0.71"^ much harder 
than the stacked signal of IRXs undetected in the 4Ms image (HR 
= -0.48 ± 0.11). This strongly suggests that the lMs CDFS 
stacked signal is indeed contaminated by hard sources which are 
now detected in the 4Ms CDFS data. Narrow band stacking of IRXs 
that are detected in the 4Ms CDFS but not the lMs CDFS (Figure(8]l 
better illustrates the full effect of the 4Ms detected population. The 
stacked spectral shape of X-ray undetected IRXs and 4Ms detected 
IRXs is similar up to 1.5 keV (approximately 4.5 keV rest-frame). 
At higher energies the stacked signals diverge and the 4Ms detected 
IRXs dominate the X-ray undetected IRXs. Including the 4Ms de- 
tected sBzKs in the stacking analysis (as in the lMs CDFS) causes 
the X-ray undetected IRXs to appear to have a much harder signal. 
We conclude that this is the main reason why th e stacked signal of 
X-ray undetected IRXs in iDaddi et al. (2007b) is harder than the 
stacked signal we find for the 4Ms CDFS X-ray undetected IRXs. 

I Alexander et alj J20 1 lh took a more in depth look at the X-ray 
properties of s BzKs, limiting their sample of sBzKs (an updated 
version of the IDaddi et alj d2007bl) sample) to those with < 5.5 
arcminute O AA and cross-c orrelating them with the 4Ms CDFS 
catalogue of IXue et alj d201ll) using a 1.5 arcsecond matching ra- 
dius. They directly detect 47 sBzKs, where 60% are IRXs and the 
remaining 40% are IRNXs. Additionally 23% (1 1/47) of the sBzKs 
are classified as heavily obscured AGN because they exhibit a hard 
X-ray spectral shape (V < 1). If we apply a similar cut in OAA to 
our sample (see Table 2) we are left with 35 X-ray detected sBzKs, 
with 51% classified as IRXs and 49% classified as IRNXs. Using 
the HR of these sBzKs we define 46% of these sBzKs as heavily 
obscured (HR > -0.08; see Section|3~TT>. 

We cross correlate the I Alexander et all d201ll) X-ray detected 
sBzKs with the 35 sBzKs with OAA < 5.5 in our sample (Table 2) 
using a search radius of 1.25 arcseconds and find that 23 sBzKs are 
common to both samples. There 19 IRXs and 14 IRNXs within this 
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Figure 8. Flux estimates for narrow band stacking of 72 X-ray undetected 
IRXs (red diamonds) and 6 IRXs that were only detected in the 4Ms CDFS 
data (black diamonds) using 4Ms CDFS images. The horizontal error bars 
represent the width of the energy band that is being stacked. 



"shared" sa mple. There are 8 sourc es which are defined as heavily 
obscured by Alexand er et alj feOllh and 9 sources which we define 
as heavily obscured in this work (see Section IXIV These heavily 
obscured sources show good agreement, with 7 of the sBzKs de- 
fined as heavily obscured in both t his work and I Alexander et al.l 
fcOllh . The X-ray fluxes from the lxile etalj d201lh 4Ms CDFS 
catalog used in the I Alexander et all d201 lh "were calculated by ex- 
tracting X-ray data over 0.5-2 keV (soft), 2-8 keV (hard) and 0.5-8 
keV (full) ranges sample and then integrating over power-laws in- 
dividually derived from the hard-to-soft band ratio (H/ S) of each 
source. We recalculated flu xes for the 23 "sh ared" sources using the 
same energy ranges as the lXueetal1 l l201ll) 4Ms CDFS catalogue 
and r = 1.4. There are inevitably some discrepancie s between 
the flu xes presented in this work and those presented bv lXue et all 
J20l 1) because of the difference in F values used to calculate flux. 
Irrespective of this the majority of fluxes show good agreement (lie 

within l a error bars). 

The Alexand er et al. I J201lh sample has 24 X-ray detected 
sBzKs which were not present in our sample, 8 of which we iden- 
tified as sBzKs but were not detected in our 4Ms CDFS catalogue, 
while the remaining 16 were not in our sBzK sample. The lXue et al.1 
J20I 1) 4Ms CDFS catalogue has 740 sources while our more strin- 
gent 4Ms catalogue has 569 sources (see Section l2"3l . so it is to be 
expected that they detect a larger number of sBzKs. The 8 galaxies 
we identify as sBzKs but do not detect in our catalogue are all low 
signifi cance detections (have low binomial detection probabilities) 
in the lXue et all ( 1201 II) 4Ms CDFS catalogue, hence likely lie just 
below the detection threshol d of our catalogue. The remaining 16 
X-ray detected sBzKs in the I Alexander et al] d201 ll) sample were 
omitted from our sBzK sample because they did not satisfy one 
or more of the criteria given in Section |2~T1 We have tried to k eep 
our sample as similar as possible to that of lDaddi et all d2007bl) but 
there are several key differences. We use a different IRAC 3.6/im 
selected catalogue to ensure the 24/im sources are deblended. Ad- 
ditionally 24/im source photo metry in our sample was extracted 
using the ConvPho t software dde Santis et al J [20071) whereas the 
iDaddi et all J2007bh sample uses a more basic fixed matching ra- 
dius method. These are the likely source of the discrepancies be- 



tween our sBzK selection. There are also 1 2 X-ray detected sBzKs 
in our sample which are not present in the I Alexander et all d201 ll) 
sample. Again this is due to the aforementioned discrepancies be- 
tween our respective X-ray point source catalogues and optical 
sam ple selection criteria. 

Alexander et al] J20T l|) s tack the remaining X -ray undetected 
sBzKs in a similar fashion to lDaddi et all d2007bh using the 4Ms 
CDFS images. They stack 47 X-ray undetected IRXs and 116 
X-ray undetected IRNXs and find effective F = 1.4 ± 0.3 and 
T = 2.0 ± 0.4 respectively. The effective F values of the stacked 
emission of X-ray undetected IRXs and IRNXs agree with those 
we calculate using our samples in 4Ms CDFS at the 1 a level (see 
Table 5). The effective F of the X-ray undete cted IRXs is also sig - 
nificantly sof ter than the valu e calcu lated in lDaddi et all d2007bl) . 
Similar to us. I Alexander et all d201 11) conclude that the majority of 
hard sources have been resolved in the 4Ms CDFS data and that 
there is no longer strong evidence that heavily obscured AGN are 
the dominant population amongst the remaining X-ray undetected 
sBzKs. 



5 DISCUSSION 

Using the deepest data available, we have investigated the X-ray 
properties of a large sample of star-forming BzK-selected galax- 
ies in the redshift range z = 1.2 - 3.0 in the CDFN and CDFS. 
The X-ray emission is particularly important when assessing the 
AGN content of these galaxies. Specifically, a direct detection 
of the X-ray emission makes the presence of an actively grow- 
ing supermassive black hole very likely, as the implied luminos- 
ity (L x > 10 42 erg s" 1 ) can only be exceeded by the most ex- 
treme star-forming galaxies. We find a relatively high (6 % ap- 
proximately), direct X-ray detection fraction for the sBzK galax- 
ies overall, indicating a very substantial AGN content. Also, con- 
sistent with previous analysis, we find a very large fraction of the 
directly-detected sBzK to be heavily obscured, based on their hard- 
ness ratios. The implied column densities are substantial (usually 
Njj > 10 23 cm -2 ). Four of the X-ray detected sBzKs have been 
identified as Compton thick AGN based on their X -ray spectra in 
previo us work: cdfs4Ms_077 and cdfs4Ms_283 by Feruglio et al. 
d201ll) and cdfs4Ms_126 and cdfs4Ms_226 by Bright man & UedJ 
d2012bl) . Of these four Compton thick AGN, 3 are IRXs and 1 is 
an IRNX. 

Stacking the emission from the sBzKs which are not detected 
directly shows a strong signal in both the soft and hard X-ray bands. 
By stacking in relatively narrow energy bands we have constructed 
a crude spectrum of these undetected objects. This is composite 
in nature, with a soft component dominating at low energies and 
harder emission at high energy. The hard spectrum above ~ 2 keV 
is a clear signature of obscured AGN activity in part of the sBzK 
population at a level below the threshold of the direct detection. 
Part of the soft emission could be scattered light from these ob- 
scured AGN. On the other hand it is also likely that unobscured, low 
luminosity unobscured AGN emission and soft X-rays from star- 
forming processes also contributes to the stack at low energies. We 
are likely therefore seeing a mixture of obscured AGN, unobscured 
AGN and star-forming galaxies. Disentangling these is rather diffi- 
cult. On the whole, however, our results support the idea that sBzK 
galaxies contain a rath er active population of obscured AGN, con - 
firming previous work dDaddi et alj2007bl [Alexander et alj|201 ll) . 

Where our results diverge somewhat from previous analysis 
(e.g. that of lDaddi et all 12007b) is when considering subsamples 
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of the sBzK population. Specifically, we have examined the direct 
detections and stacked X-rays from subsamples based on their "IR 
Excess", the ratio of the total star formation rate to that of the dust- 
corrected UV SFR. IRXs do exhibit a significantly higher X-ray 
detection fraction (by a factor 3 — 4) than IRNXs. However, the 
majority of this can be relatively simply explained given the ex- 
tremely strong correlation between the IR excess parameter and 
the total infrared luminosity. When controlling for IR luminosity, 
there is no statistically-significant difference in the X-ray detected 
fraction comparing IRXs and IRNXs. 

Stacking results for the fainter sources confirm this effect. The 
stacked IRXs give a slightly higher flux, but the hardness ratios of 
the IRXs and IRNXs are very similar. Stacking in narrow bands 
also shows that the two populations have very similar spectra to 
each other. In particular, both stacks show a hard component in- 
dicating the presence of obscured AGN in some of the objects, 
too faint to be detected individually. Overall, it seems clear that 
the IR excess parameter is not decisive in selecting obscured AGN 
or AGN more generally from the sBzK population. The fact that 
the IR-luminous sBzKs show a high X-ray detection fraction could 
be attributable to a number of factors. First the BzK selection will 
tend to yield samples of galaxies with high stellar mass and red 
colour s, which are known to preferential ly host luminous AGN 
(e.g. lNandra et al.ll2007l : lBundv et alj2008t) . Galaxies with high IR 
luminosities will also tend to be the most luminous and massive 
sources themselves. The high IR luminosities and star formation 
rates also suggest copious amounts of gas and dust, which could 
account for the high levels of obscuration in these sources. On the 
other hand, it is puzzling why the presence of an AGN as inferred 
by the presence of (hard) X-ray emission apparently bears so lit- 
tle relation to excess 24/^m emission. Apparently there are many 
galaxies without a strong hot dust component which nonetheless 
must host an AGN based on their X-ray emission and, similarly, 
many which have very strong 24/j.m but without a powerful AGN. 

iDaddi et"al] ( l2007bT) have suggested that the obscuration levels 
in the IRXs may in fact be so high that they are "Compton thick" 
AGN, where even the X-ray emission is suppressed and observed 
only in scattered light. If so a considerable amount of additional 
black hole growth may be hidin g in these objects. This issue has 
been investigated further by Alexander et al. 20 lj] who examined 
the X-ray spectral properties of directly detected sBzKs using the 
4Ms CDFS data. Their analysis resulted in a sharp downward revi- 
sion in t he estimated Comp ton thick AGN fraction in BzKs com- 
pared to lDaddi et all |2007b), who used X-ray stacking analysis to 
show that the hard X-ray emission of the IRXs resembled the typi- 
cal spectrum of Compton thick AGN, with a very hard spectrum at 
high energies. 

Our stacking results confirm the presence of a hard high en- 
ergy spectrum in the sBzK population. Overall the spectrum of the 
stacked sources is of a "composite" nature, with soft emission ei- 
ther from star formation or unobscured AGN, and a hard signal 
from obscured AGN. Determining whether the obscuration in the 
obscured AGN is Compton thin or Compton thick is extremely dif- 
ficult, as the X-ray continuum from moderately obscured Comp- 
ton thin AGN (A?h ~ 10 23 cm -2 ) and the reflection dominated 
spectrum of Compton thick AGN is quite similar. The telltale sig- 
nature is the intense iron Ka line (E = 6.4keV rest-frame) pre- 
dicted in the Compton thick case. This signature has been found 
in at least one sBzK gala xy dFeruglio et "al] |201lh and is inferred 
by lAlexander et al.l J20 1 lh to be present in a handful of additional 
sBzKs, based on the co-addition of the X-ray spectra of several di- 
rectly detected sBzKs with hard X-ray colours in the 4Ms CDFS 



data. As discussed ab ove the conclusion o f widespread Compton 
thick AGN in BzKs bv lDaddi et ail (12007b) was based on the shape 
of the stacked spectrum. With a much larger sample of BzK galax- 
ies, and far deeper X-ray data, it may have been hoped that our 
stacking analysis would reveal the iron Ka line predicted in the 
Compton thick AGN case. In the event, however, our results nei- 
ther confirm not deny the presence of intense iron Ka emission 
in the galaxies. The primary issue is that the deeper X-ray data 
are able to resolve out the brighter X-ray sources from the stack, 
which carry the majori ty of the signal-to-noise ratio. As noted by 
lAlexander et alj ( l201ll) . and confirmed by us, a great number of 
the objects directly detected in the 4Ms CDFS C handra data, but 
not in the lMs data used bv lDaddi et ail d2007bl) , have very hard 
spectra. These sources are undoubtedly obscured, and a fraction 
of them may be Compton thick, but once removed from the stack 
evidence for a Compton thick AGN in the residual sBzK popula- 
tion is weakened, rather than strengthened. We conclude that it is 
particularly difficult to base estimates of the contribution of Comp- 
ton thick AGN on a stacking analysis, as stacking will result in a 
mixture of X-rays from obscured AGN, unobscured AGN and stel- 
lar processes, and can be strongly affected or even dominated by 
a relatively small number of sub-threshold sources. Definitive con- 
clusions regarding the presence of the tell-tale iron Ka emission 
in faint Compton thick AGN awaits future experiments with much 
larger collecting are a than Chandra, such as the Athena concept 
dBarcons et all2012l) . 

In addition to impacting on our knowledge of the AGN content 
of sBzK galaxies, our X-ray results also have interesting implica- 
tions about the nature of the IR-excess sBzK population. The basis 
of the X-ray analysis of the IRXs is that the additional emission at 
24/im could be due to hot dust emission powered by a buried AGN. 
Our analysis has clearly shown, however, that the AGN content 
(and obscuration properties) of IRXs and IRNXs are very similar. 
This raises the question of the origin of the large discrepancy be- 
tween the 24/^m-derived star formation rates and the dust-corrected 
UV measurements. As these are apparently not due to AGN activity 
it indicates that the derivation of one or other of the SFRs using the 
methodology adopted in the current paper is faulty. The discrep- 
ancy is largest at high IR luminosity, where the 24/^m SFR greatly 
exceeds the UV estimates, in extreme cases by 2 orders of magni- 
tude. This implies either that the 24/im SFR is overestimated, the 
UV underestimated, or both. Recent results from Herschel suggest 
that both effects could be in play. Firstly, 24/im fluxes appear sys- 
tematically Jo_owr;esthn^e the_SF^^ 1 £sp£cMly of high redshift 
galaxies dElbaz et alj|201 ll : iNordon et alj|2012h . Secondly, at the 
highest IR luminosities and thus SFRs, the lines of sight to the star 
forming regions may become complet ely opaque to UV em ission, 
being obscured by optically thick dust JGoldader et alj .2002). This 
makes the rest-frame UV emission, even if corrected for dust ex- 
tinction, an extremely unreliable estimate of the star formation rate 
and specifically it will lead one to underestimate the SFR. 

A final interesting point from the current work stems from the 
X-ray analysis of the passive BzK galaxies. While our sample of 
pBzK is not large (26 objects) it is sufficient to be able to make the 
first estimate of the AGN content of these galaxies. Of 18 pBzKs 
in the redshift range z — 1.2 — 3.0 we find two X-ray detections. 
Expanding out the full redshift range spanned by our BzK sample 
the detections rise to 4/26 objects. While the number statistics are 
still small, the X-ray detection rate and AGN content of the pBzKs 
is at least consistent with and perhaps even slightly higher than 
detection rate in sBzKs. The high AGN fraction in sBzKs has been 
interpreted bv lDaddi et aIH 2007b) as evidence for rapid coeval star 
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formation and black hole growth in this galaxy class. It appears, 
however, that this is not robust, as BzK galaxies without rapid star 
formation also show an unusually high fraction of AGN. So while 
a large amount of black hole growth in galaxies at z ~ 2 does 
indeed occur in galaxies that are also actively growing their stellar 
component, this is arguably more likely to be due to the fact that 
most massive galaxies are high z ~ 2 are star forming galaxies. 
It does not necessarily present a causal and or astrophysical link 
between black hole growth and star formation. 



6 SUMMARY 

Using 4Ms CDFS and 2Ms CDFN Chandra/ ACIS imaging we have 
investigated the X-ray properties of BzK-selected galaxies at z ~ 
2. The AGN content and obscuration properties have been studied 
using a combination of direct detection, broad band stacking and 
narrow band stacking analyses. The main results are listed below: 

(i) We directly detect 5.6 ± 0.6% (81/1445) of star-forming 
sBzKs in X-rays, finding 49 sBzKs in the 4Ms CDFS and 32 sBzKs 
in the 2Ms CDFN Chandra ACIS-I images. The direct detection 
limits are such that almost all of these will host AGN. We di- 
rectly detect a greater proportion of IR Excess sBzKs relative to 
IR Non-Excess sBzKs. However, this can be attributed entirely to 
the fact that the IR excess criterion also selects the most IR lumi- 
nous sBzKs. The X-ray detected fractions of IR Excess sBzKs and 
IR Non-Excess sBzKs are very similar when binned by Ljr, 

(ii) A high proportion (44%) of X-ray detected sBzKs are heav- 
ily obscured AGN (HR > —0.08), while the remaining sBzKs 
are relatively unobscured AGN, with a few perhaps purely star- 
forming. The IR Excess sBzK and IR Non-Excess sBzK subsam- 
ples exhibit very similar proportions of heavily obscured AGN. 
Four of the X-ray detected sBzKs have been identifi ed as Compton 
thick A GN from their X-ray spectra in the work of iFeruglio et al.l 
d201 lh and lBrightman & Uedal l l2012bh . 

(iii) Broad band stacking of 688 X-ray undetected sBzKs, as 
well as IR Excess sBzK and IR Non-Excess sBzK subsamples, in 
the 4Ms CDFS and 2Ms CDFN yields soft stacked signals with 
effective F ~ 2. Stacking of sBzKs binned by redshift and Lir 
produces similar results. 

(iv) Narrow band stacking of sBzKs reveals a "composite" spec- 
trum with soft emission at low energies and hard emission at high 
energies. It is not possible to disentangle the emission from heav- 
ily obscured AGN, unobscured AGN and star-formation within the 
stacked signal. The resolution of the narrow band stacking is not 
sufficient to confirm or deny the presence of Iron Ka line emis- 
sion, so it cannot be established definitively whether or not a sig- 
nificant population of Compton thick AGN reside in the sBzKs. 
The IR Excess sBzKs and IR Non-Excess sBzKs have almost iden- 
tical spectral shapes and, as with broad band stacking, there is no 
evidence that Compton thick AGN are more prevalent within IR 
Excess sBzKs than IR Non-Excess sBzKs. 

(v) The stacked emission of IR Excess sBzKs in the lMs CDFS 
is dominated by 6 IR Excess sBzKs lying just below the detection 
threshold at this depth, and these were subsequently detected in the 
deeper 4Ms CDFS. These sources are on average much harder (ef- 
fective T = 0.7±°:?) *an the X -ray undetected sBz Ks and are 
respons ible for the hard stacked signal observed by Daddi et al. 
d2007bh . 

(vi) Our results show that on the whole the discrepancy between 
24^m and UV-derived star formation in the BzK population is not 
due to AGN activity enhancing the 24/im flux. At high Lir (hence 



SFR) the 24^im overestimates the SFR, and the UV underestimates 
SFR for these BzK galaxies. The apparent lack of correspondence 
between X-ray emission from an AGN and 24/im excess is a puz- 
zle which needs to be investigated in detail on an object-by-object 
basis. 

(vii) Passively evolving BzKs have an even higher X-ray detec- 
tion rate (~ 15%) than sBzKs (~ 7%) and hence a high fraction 
of AGN. This is contradictory to the idea that black hole growth at 
z ~ 2 is coeval with star-formation. The connection between the 
two appears to be purely circumstantial, rather than causal. 
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Table 5. Broad band stacking of X-ray undetected sBzK and pBzK galaxies in 4Ms CDFS and 2Ms CDFN. All fluxes calculated using flux conversion assuming T = 1.4. Column(l): Sample of galaxies stacked; 
column (2): No. of galaxies used stacked; column (3): soft band S/N; column (4): hard band S/N; column (5): full band S/N; column (6): 0.5-2keV flux, derived from soft band stacked counts, units 10 -17 erg cm~ 2 
s _1 , if we assume T = 2.0 the soft band fluxes would be systematically 10% lower; column (7): 2-10keV flux, derived from hard band stacked counts, units 10~ 17 erg cm~ 2 s — 1 ; column (8): 0.5-10keV flux, 
derived from full band stacked counts, units 10~ 17 erg cm~ 2 s _1 ; column (9): HR, HR= (H — S)/(H + S), where H and S are hard- and soft-band count rates respectively; column (10): effective F of stack, 
derived from HR. 



Sample 


No. Sources 




S/N 2 _7fceV 


S/No.5-7fceV 


Fo.5-2keV 


^2-lOfceV 


Fl).5-10keV 


HR 






r 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 






(10) 


4Ms CDFS sBzKs 


370 


16.00 


4.24 


12.97 


0.51 ± 0.03 


0.90 ±0.21 


1.71 ± 0.13 


-0.53 ± 0. 


08 


2. 


1 +0.3 
— 0.2 


4Ms CDFS IRXs 


72 


10.64 


3.42 


9.33 


0.82 ± 0.08 


1.66 ±0.49 


2.90 ± 0.31 


-0.48 ± 0. 


11 


2 


q+0.3 
u -0.3 


4Ms CDFS IRNXs 


298 


12.40 


3.02 


9.75 


0.43 ± 0.03 


0.71 ±0.24 


1.43 ± 0.15 


-0.55 ± 0. 


11 


2. 


n+0.4 
-0.3 


2Ms CDFN sBzKs 


318 


12.88 


4.17 


11.50 


0.59 ± 0.05 


1.22 ±0.29 


2.10 ± 0.18 


-0.47 ±0. 


09 


2. 


u -0.2 


2Ms CDFN IRXs 


107 


10.34 


3.31 


9.34 


0.88 ± 0.09 


1.69 ±0.51 


3.08 ± 0.33 


-0.50 ±0. 


12 


2. 


n +0.4 
u -0.2 


2Ms CDFN IRNXs 


211 


8.20 


2.74 


7.34 


0.44 ± 0.05 


0.98 ±0.36 


1.61 ± 0.22 


-0.45 ± 0. 


15 


1 


q+0.5 
9 -0.3 


Combined sBzKs 


688 


20.44 


5.77 


17.07 


0.55 ± 0.02 


1.04 ±0.18 


1.89 ± 0.11 


-0.50 ±0. 


06 


2. 


n +0.2 
u -0.1 


Combined IRXs 


179 


14.81 


4.73 


13.12 


0.86 ± 0.06 


1.68 ±0.36 


3.00 ± 0.23 


-0.49 ± 0. 


08 


2 


n+0.3 
u -0.2 


Combined IRNXs 


509 


14.84 


3.95 


12.08 


0.43 ± 0.03 


0.82 ±0.20 


1.50 ±0.13 


-0.51 ± 0. 


09 


2. 


1+0.4 
-■--0.2 


Combined pBzKs 


13 


2.99 


2.09 


3.47 


0.49 ± 0.16 




2.23 ± 0.62 









Table 6. Broad band stacking of redshift and IR luminosity binned X-ray undetected sBzKs. Galaxies in 4Ms CDFS and 2Ms CDFN have been combined to boost the S/N ratio of the stacks. All fluxes calculated 
using flux conversion assuming T ~ 1.4. Column(l): Sample of galaxies stacked; column (2): No. of galaxies used stacked; column (3): soft band S/N; column (4): hard band S/N; column (5): full band S/N; column 
(6): 0.5-2keV flux, derived from soft band stacked counts, units 10~ 17 erg cm~ 2 s — 1 ; column (7): 2-10keV flux, derived from hard band stacked counts, units 10~ 17 erg cm~ 2 s~ 1 ; column (8): 0.5-10keV flux, 
derived from full band stacked counts, units 10~ 17 erg cm~ 2 s _1 ; column (9): HR, HR= (H — S)/(H + S), where H and S are hard- and soft-band count rates respectively; column (10): effective F of stack, 
derived from HR. 



Sample No. Sources S/N . 5 _ 2 fceV S/N 2 _ 7 fe e v S/N . 5 -7fceV Fo.5-2keV ^2-lOfceV ^o.S-lOfceV HR T 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



1.2 < z < 1.8 


340 


12.62 


4.06 


10.91 


0.38 ± 0.03 


0.71 ±0.18 


1.29 ±0.12 


-0.47 ±0.10 


2 0+ ' 3 


1.8 < z < 2.4 


304 


14.28 


4.12 


11.97 


0.47 ±0.03 


0.84 ±0.19 


1.54 ±0.13 


-0.48 ±0.09 


2 +0 ' 3 


2.4 < z < 3.0 


97 


6.46 


0.62 


4.31 


0.36 ±0.05 




1.01 ±0.21 






Li a < 10 11 


108 


5.96 


2.12 


5.31 


0.33 ±0.06 


0.71 ±0.33 


1.16 ±0.21 


-0.42 ±0.20 


1 S+ ' 7 


10 11 < L IR < 10 12 


393 


11.16 


3.60 


9.46 


0.29 ±0.03 


0.64 ±0.16 


1.01 ±0.11 


-0.40 ±0.11 




10 12 < L IR < 10 13 


203 


13.42 


3.73 


11.31 


0.56 ±0.04 


0.86 ±0.23 


1.78 ±0.16 


-0.54 ±0.10 


2 1+ ' 4 


Lir > 10 13 


59 


13.44 


7.03 


14.14 


1.21 ±0.09 


2.89 ±0.47 


4.40 ±0.33 


-0.37 ±0.08 


i 7 +o.2 
1 -'-o.i 



Table 7. Broad band stacking of X-ray undetected sBzK galaxies in lMs CDFS and sBzKs that were subsequently detected in the 4Ms CDFS. Galaxies in 4Ms CDFS and 2Ms CDFN have been combined to boost 
the S/N ratio of the stacks. All fluxes calculated using flux conversion assuming T = 1.4. The 4Ms detected galaxies i*b.5-i0fceV is lower than Fi-WkzV because the powerlaw (r=1.4) used to approximate the 
spectra is too soft. Column(l): Sample of galaxies stacked; column (2): No. of galaxies used stacked; column (3): soft band S/N; column (4): hard band S/N; column (5): full band S/N; column (6): ().5-2keV flux, 
derived from soft band stacked counts, units 10~ 17 erg cm~ 2 s _1 ; column (7): 2-l()keV flux, derived from hard band stacked counts, units 10~ 17 erg cm~ 2 s _1 ; column (8): 0.5-10keV flux, derived from full 
band stacked counts, units 10 -17 erg cm -2 s~ x ; column (9): HR, HR= (H — S) / (H + S), where H and S are hard- and soft-band count rates respectively; column (10): effective T of stack, derived from HR. 



Sample 


No. Sources 


S/No.5_2fceV 


S/N2-7feeV 


S/N .5_7fe e v 


Fo.5-2keV 


F'2-10keV 


Fo.5-10keV 


HR 


r 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


lMs CDFS sBzKs 


391 


9.32 


3.64 


8.84 


0.51 ±0.06 


1.29 ±0.36 


1.94 ±0.22 


-0.39 ±0.12 


1 S+ ' 3 


lMs CDFS IRXs 


80 


5.56 


3.60 


6.32 


0.70 ±0.13 


2.96 ±0.82 


3.20 ±0.51 


-0.15 ±0.16 




lMs CDFS IRNXs 


311 


7.58 


2.17 


6.62 


0.46 ±0.06 


0.86 ± 0.40 


1.62 ±0.24 


-0.51 ±0.17 


9 1+O.6 


4Ms CDFS detected sBzKs 


12 


9.69 


9.96 


13.74 


2.41 ±0.25 


16.02 ± 1.60 


13.89 ± 1.01 


0.07 ±0.07 


0.8l° ; 2 


4Ms CDFS detected IRXs 


6 


7.52 


8.57 


11.29 


2.72 ±0.36 


20.61 ±2.40 


16.98 ± 1.50 


0.14 ±0.09 


0.7t°; 2 


4Ms CDFS detected IRNXs 


6 


6.13 


5.27 


7.37 


2.11 ±0.35 


11.15 ±2.12 


10.80 ± 1.36 


-0.04 ±0.12 


i-ol8:I 



